Crossover behavior in a mixed-mode fiber bundle model 
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We introduce a mixed-mode load sharing scheme in fiber bundle model. This model reduces exactly 
to equal load sharing (ELS) and local load sharing (LLS) models at the two extreme limits of a 
single load sharing parameter. We identify two distinct regimes: a) Mean-field regime where ELS 
mode dominates and b) short range regime dominated by LLS mode. The crossover behavior is 
explored through the numerical study of strength variation, the avalanche statistics, susceptibility 
and relaxation time variations, the correlations among the broken fibers and their cluster analysis. 
Analyzing the moments of the cluster size distributions we locate the crossover point of these regimes. 
We thus conclude that even in one dimension, fiber bundle model shows crossover behavior from 
mean-field to short range interactions. 



I. INTRODUCTION 

Fiber bundle model represents a simple, 
stochastic fracture-failure process Q in mate- 
rials subjected to external load. The model 
consists of three basic ingredients: (a) a 
discrete set of N elements located at sites 
of a lattice (b) a probability distribution of 
the strength threshold of individual elements 
(fibers) (c) a load-transfer rule which dis- 
tributes the terminal load carried by the failed 
fibers to the surviving fibers. The model study 
was initiated by Peirce in the context of 
testing the strength of cotton yarns. Since 
then this model has been studied and modi- 
fied by many groups [3-25] using analytic as 
well as numerical methods. Fiber bundles are 
of two classes with respect to the time depen- 
dence of fiber strength threshold: 'Static' bun- 
dles contain fibers whose threshold strengths 
are independent of time and such bundles are 
subjected to quasi static loading, i.e., the load 
is increased steadily up to the complete fail- 
ure of the bundles. The load or stress a (load 
per fiber) is an independent variable here and 
the strength of the bundle is determined by 
the maximum value of the applied load or 
stress (<7 C ) that can be supported by the bun- 
dle. On the other hand 'dynamic' bundles con- 
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sist of fibers having time dependent strength 
and the fibers fail due to fatigue 0, 0, @, u} 
after a period of time which varies fiber to 
fiber. The time taken for complete failure 
is called the lifetime of the bundle. Accord- 
ing to the load sharing rule, fiber bundles are 
being classified into two groups: Equal load- 
sharing (ELS) bundles [8-17] or democratic 
bundles and local load-sharing (LLS) bundles 
E 13 ED| In the ELS models all the in- 
tact fibers equally share the terminal load of a 
failed fiber, whereas in LLS model the termi- 
nal load gets shared among the intact nearest 
neighbors. ELS models show phase transition 
from partial failure to total failure at a critical 
strength (cr c ). The critical behavior in the fail- 
ure dynamics of ELS bundles has been solved 
analytically 0,0] and the universality of the 
ELS model has been established recently. 
However the strength of LLS models goes to 
zero 0, |U E3| at the limit of infinite sys- 
tem size and this does not permit any critical 
behavior in the failure process. 

The ELS and LLS models belong to two 
opposite extremes with respect to the spatial 
correlations in stress redistributions. These 
models do not incorporate any type of stress 
gradient among the intact fibers which is an 
usual expectation. Therefore a load sharing 
scheme in between ELS and LLS should be a 
realistic approach to study the failure of het- 
erogeneous materials. Hansen and Hemmer 
[2J| introduced a 'A model' to interpolate be- 
tween ELS and LLS models where A is an ad- 
justable stress transfer factor. Although they 
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conjectured the existence of a critical crossover 
value A c which separates the mean field (ELS) 
regime and the short range (LLS) regime, what 
would be the exact crossover point was not an- 
swered. A recent approach by Hidalgo et al 
(2^ | incorporate both the ELS and LLS mode 
introducing an effective range of interaction 
parameter (7) which is actually the power of 
the stress redistribution function. They ob- 
served crossover behavior in strength variation 
and in the avalanche statistics of the failures. 
Also they determined the crossover point (7 C ) 
through the moment analysis of the cluster size 
distributions before total failure. 

In this paper we develop a mixed-mode load 
sharing (MMLS) model which interpolates the 
ELS and LLS models correctly. We intend to 
study whether this model shows a continuous 
transition from mean-field (ELS) behavior to 
extreme statistics (LLS), or there exists a def- 
inite crossover point. 

We organize this paper as follows: After in- 
troduction (section I) we present our MMLS 
model in section II. Section III contains the ob- 
servations of Crossover behavior through nu- 
merical study of the model. The analysis to 
determine exact crossover point is given in sec- 
tion IV. The final section (V) is devoted for 
discussions including our conclusions. 



II. THE MODEL 

Our mixed-mode load sharing (MMLS) 
scheme is basically a coupling of ELS and LLS 
mode: When a fiber fails, a fraction (g) of 
its terminal load gets shared among the near- 
est neighbors of the failed fiber (LLS rule) 
and the rest (1 — g fraction) is distributed 
equally among all the surviving fibers (ELS 
rule). Here 'g' is the weight parameter of the 
MMLS scheme. Therefore, the model reduces 
exactly to ELS model for g = and for 5=1, 
it becomes pure LLS one. As we have cho- 
sen 1 — D fiber bundle model (with periodic 
boundary condition), the number of nearest 
neighbors is always two. We study the behav- 
ior of the model for the entire range < g < 1 
using Monte-Carlo simulations for step-wise 
equal load increment (l^Lll^.lr^ | until the total 
failure of the bundle. During the entire study 



we consider uniform (on average) distribution 
of fiber strength threshold in the bundle. 



III. THE CROSSOVER BEHAVIOR 

A. Strength of the bundle 

It is known since Daniels 0] that the ELS 
bundles have a nonzero strength (cr c ) above 
which the bundle fails completely. Recently 
it has been shown analytically [lfiL Ilf3 | that 
for uniform fiber threshold distribution, bun- 
dle's strength approaches the value 1/4 as sys- 
tem size goes to infinity. On the other hand 
LLS bundles do not have any nonzero strength 
HHH El- In our MMLS model we intend to 
study the strength variation of bundles with 
system size as well as with the weight param- 
eter g. 



0.26 
0.24 
0.22 
0.2 

D° 0.18 
0.16 
0.14 
0.12 

0.1 



N=5000 
N=10000 
N=20000 
N=50000 
N=1 00000 



1.2 1.4 



0.2 0.4 0.6 0.8 

g 

Fig. 1: The strength of the bundle for differ- 
ent system sizes (N) as a function of the weight 
parameter g. 

As g increases, the bundle becomes weaker 
due to the short-range (LLS) interactions. 
Therefore <j c decreases with increasing g val- 
ues (Fig. 1). We can see that a c seems to 
be independent on system size (dominance of 
ELS) up to g = 0.7 and beyond g — 0.8, a 
strong system size dependence (dominance of 
LLS) appears. This observation is supported 
by Fig. 2, where we have shown the logarith- 
mic size dependence of a c . Up to g = 0.7, the 
curves eventually become flat as the system 
size increases. But for g > 0.8, all the curves 
fall (following inclined straight line) . Thus the 



3 



two regimes are differentiated clearly. 
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Fig. 2: The logarithmic size dependence of 
bundle's strength for different values of g. 
straight lines represent the best fit. 
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20000). The dotted line represents mean-field re- 
sult having exponent value —5/2. Clearly, the up- 
per group of curves can be fitted by the mean-field 
power law whereas the lower group does not show 
power law at all. 



Here we have measured (Fig. 3) the 
avalanche size distributions for different g val- 
ues. Clearly, two groups of curves appear. 
The upper group (0 < g < 0.7) can be fit- 
ted with the mean-field result (—5/2) where 
as the lower group (0.8 < g < 1.0) show a 
clear deviation from the power law. 



B. Avalanche size distribution 

The avalanche size distribution character- 
izes the fracture process by reflecting the 
precursory activities toward complete failure. 
This can be related to the acoustic emissions 
observed in material failure \'2(i 1271 l28j . Hem- 
mer and Hansen showed @ analytically that 
for ELS models the avalanche size distribu- 
tion follows an universal power law with expo- 
nent value —5/2. But for LLS models the nu- 
merically estimated apparent exponent value 
is quite larger 4.5 Later it has been shown 
analytically (for flat threshold distributions) 
that for LLS model, no universal power-law 
asymptotics exists [HI ]. 
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Fig. 3: Avalanche size distribution for dif- 
ferent values of the weight parameter g (averag- 
ing over 5000 configurations for system size N = 



C. The susceptibility and relaxation time 
variations 



Recently, the dynamic response parameters, 
susceptibility (x) Hp. H3. \wi flfl l2fl| and re- 
laxation time (t) \WL iW^ have been studied in 
fiber bundle models. The susceptibility is de- 
fined as the number of fibers fail due to an in- 
finitesimal change of the external stress (a) on 
the bundle and the relaxation time is the time 
(number of stress redistributions) the bundle 
takes to come to a stable fixed point at an ex- 
ternal stress (a). For ELS model, the suscepti- 
bility and relaxation time seem to follow power 
law with the applied stress and both of them 
diverge [lj, Ha U3 at the critical strength 
a c : x ~ (c c - <t)" 1/2 and r ~ (a c - ct) _1/2 . 
However, one can not expect such scaling be- 
havior in LLS models due to the absence of 
'critical' strength. The step-wise equal load 
increment method [TH enables to measure 
X and t for different values of g (Fig. 4). The 
power law behavior (with mean-field exponent 
— 1/2) remains unchanged up to g = 0.7 and 
for g > 0.8 the curves do not follow power laws 
at all. Thus the susceptibility and relaxation 
time variations also suggest a transition from 
the mean-field to short range behavior to hap- 
pen in between g = 0.7 and g = 0.8. 
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(CT-C) (O c -0) 
Fig. 4: The susceptibility (x) and relaxation 
time (r) variations for different g values. The 
bundle contains 10000 fibers and the data are 
averaged over 10000 configurations. 



D. Correlations among the broken fibers 




Fig. 5: The space-time diagram of the break- 
down sequence in MMLS model. The positions of 
the fibers are marked on the x axis and y axis is 
a 'time' axis where time indicates the number of 
stress redistribution starting from initial loading. 
The white color represents intact fibers while the 
black regions stands for the broken fibers. 



The breakdown sequence reflects the corre- 
lations of the breaking process While 
the ELS model simply ignores the spatial ar- 
rangement of the fibers, LLS model gives much 
importance on it. Therefore, as g increases 
(LLS mode dominates) the breaking process 
becomes more and more correlated (Fig. 5). 
Here also we can identify two distinct regimes. 
We cannot see any spatial correlation among 
the broken fibers (except near the total fail- 
ure) up to g = 0.7, whereas for g > 0.8 strong 
correlations (black patch) develop long before 
the total failure. 
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IV. DETERMINATION OF THE 

EXACT CROSSOVER POINT 
THROUGH CLUSTER MOMENT 
ANALYSIS 



The fracture process can also be character- 
ized by analyzing the clusters of broken fibers 
just before complete failure [l^. l2fl|. The size 
distributions of the clusters (n(s) vs. s) are 
shown (Fig. 6) for different values of g. Al- 
though the distributions appear as two groups, 
it is not possible to identify the exact crossover 
oint from this. Therefore we go for the mo- 
ment analysis: the A:-th moment of the cluster 
.istributions is defined 1 25 1 as 



nik = / s k n(s)ds 



(1) 



position 



Clearly the zero-th moment (mo) gives the 
total number clusters and the first moment 
gives the total number of broken fibers. We 
can get the average cluster size dividing the 
second moment (m.2) by first moment (mi). 
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tioe. As g increases the stress redistribution 
becomes more and more localized in the neigh- 
borhood of the failed fibers and a few isolated 
crack can trigger the complete rupture through 
growth and coalescence mechanism. There- 
fore the pure ELS mode can store the maxi- 
mum crack (cluster) and this capacity should 
decrease with the increase of g. We can see 
(Fig. 7) that both mo and mi decrease with 
increasing g value and they fall drastically in 
between g = 0.7 and g — 0.8. This crossover 
is much robust in case of m.2 and the average 
cluster size ^/mi), both of which show a 
: arp peak, which indicates the dominance of 
LLS mode over the ELS mode 25]. 



Fig. 6: Cluster size distributions of broken 
fibers (just before complete failure) for different 
g values (averaging over 5000 samples for N = 
20000). 




Fig. 7: The moments of the cluster size dis- 
tributions as a function of the weight parameter g 
(averaging over 5000 samples for N = 20000) . 

In case of pure ELS mode (g — 0), we 
have only long-range interaction and the clus- 
ters are randomly distributed within the lat- 
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Fig. 8: System size dependence of the 
crossover point. 

To check how the crossover point changes 
its position with system size, we have done 
the similar cluster moment analysis for sev- 
eral system sizes. We observe a weak system 
size dependence of the peak position i.e., the 
crossover point (Fig. 8). With proper extrap- 
olation we determine the crossover point (g c ) 
to be at g = 0.79 ± 0.01 for the infinite system 
size. 



V. CONCLUSION 

The fracture and breakdown of loaded ma- 
terials is basically a cooperative phenomenon 
guided by the load redistribution mechanism. 
Here 'crack' opens up when an element (fiber) 
fails after external loading. This single fiber 
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failure should affect the neighbors much than 
the distant elements (like in electric Fuse Mod- 
els Therefore a high stress concentra- 
tion (after the load redistribution) around a 
crack (failed fiber) is a natural expectation. 
The ELS models do not incorporate any spa- 
tial correlations and exhibit perfect democracy 
(mean-field), whereas the LLS models confine 
themselves within the nearest neighbor inter- 
actions. In this situation attempts 0, H3| 
to study the failure behavior in between ELS 
and LLS regimes, would be most welcome. 
Also, a recent experiment on loaded wood- 
fiber 30j demands an intermediate load shar- 
ing scheme to explain the observed strength 
variation. The 'A model' |24| becomes a LLS 
model at A = 1. But it cannot be reduced to 
a pure ELS model at A = 0, as the neighbors 
of the just broken fibers become 'immunized' 
against failure. Alth oug h the 'variable range 
of interaction' model (23 determines the exact 
crossover point, it remains silent about the sys- 
tem size dependence of crossover point, which 
is nevertheless an important issue. 

Our mixed-mode load sharing (MMLS) 
model exactly reduces to ELS model at g = 
and to LLS model at g = 1. We estab- 
lish numerically that the MMLS model in one 
dimension shows a distinct crossover behav- 
ior from mean-field to short-range interaction. 
The strength (<7 C ) variation of the bundle with 



system size, the avalanche statistics and the 
failure dynamics (susceptibility and relaxation 
time) suggest that the crossover point (g c ) 
must be in between g = 0.7 and g = 0.8. 
The cluster size analysis determines the exact 
crossover point in one dimension for several 
system sizes and a proper extrapolation sug- 
gests the crossover point to be g c = 0.79 ±0.01 
at the limit of infinite system size. For g < 
g c the model exhibits critical behavior (sup- 
ported by the power laws) for the dominance of 
ELS mode. But the fluctuations suppress any 
critical behavior after g — g c , where extreme 
statistics [lj dominates. We should mention 
that as the ultimate strength (cr c ) of the bun- 
dle continuously decreases with the increas- 
ing g value, we cannot exclude the possibil- 
ity of different critical behavior for g = 0, 
< g < g c and g > g c in higher dimensions, 
like in case of 2 — D Ising systems with dis- 
order [^lj ■ Therefore we expect this crossover 
behavior in MMLS model to be more promi- 
nent in higher dimensions. 
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